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The spin precession in a cylindrical semiconductor nanowire due to Rashba spin-orbit coupling has 
been investigated theoretically using an InAs nanowire containing a surface two-dimensional electron 
gas as a model. The eigenstates, energy-momentum dispersion, and the energy-magnetic field dis- 
persion relation are determined by solving the Schrodinger equation in a cylindrical symmetry. The 
combination of states with the same total angular momentum but opposite spin orientation results 
in a periodic modulation of the axial spin component along the axis of the wire. Spin-precession 
about the wires axis is achieved by interference of two states with different total angular momentum. 
Because a superposition state with exact opposite spin precession exists at zero magnetic field, an 
oscillation of the spin orientation can be obtained. If an axially oriented magnetic field is applied, 
the spin gains an additional precessing component. 



I. INTRODUCTION 

Semiconductor nanowires are almost ideal objects 
for studying quantum effects and electron interference 
phenomena. The use of the bottom-up approach for 
nanowire growth simplifies the preparation substantially 
and allows us to create novel confinement schemes, such 
as axial or radial heterostructuresP^ The large surface- 
to-volume ratio of nanowires means that surface prop- 
erties are crucial for discussions of transport properties, 
so that low band-gap semiconductors, e.g. InAs, InN, 
or InSb, are particularly interesting. In these systems, 
the Fermi level at the surface is pinned inside the con- 
duction bandP and an accumulation layer is formed. 
This guarantees that the conductance is sufficiently large 
even at low nanowire radius. The presence of the sur- 
face accumulation layer means that a tubular conduct- 
ing channel is formed, and this shape of the conductor 
has important implications for the magnetoconductance 
of the nanowires. An example is the theoretical predic- 
tion and experimental confirmation of flux-periodic oscil- 
lations in nanowires with a magnetic field applied along 
the wire axis The electronic states of a cylindrical two- 
dimensional electron gas in a transverse magnetic field 
were calculated by Ferrari et a/., 6 7 while Magarill et alW^ 
discussed the kinetics of electrons in a tubular conductor. 

Many concepts have been developed for planar semi- 
conductor layer systems that make use of the spin de- 
gree of freedom for device structures. The best-known 
example is the spin field-effect transistor ESHUI w hich 
uses the gate-controlled spin-precession induced by the 
Rashba effect EHUD The Rashba spin-orbit coupling orig- 
inates from a macroscopic electric field in an asymmet- 
ric quantum well. 16 Meanwhile, research activities have 
been extended to planar quasi one-dimensional struc- 
tures, which promise a superior spin controlJ^^ The 



energy spectrum and spin precession in these structures 
are governed by the interplay between confinement and 
energy splitting due to spin-orbit coupling! 20 * 21 ! Only a 
few theoretical investigations have dealt with the effect 
of spin-orbit coupling in cylindrical conductors on the 
electronic states and on the quantum transport! 4 * 9 * 22 * 23 ^ 
The spin-dynamics in curved two-dimensional electron 
gases was discussed by Trushin and Schliemann 24 while 
the weak antilocalization effect in cylindrical wires was 
studied by Wenk and KettemannP^The presence of spin- 
orbit coupling was confirmed for In]>P^and InAs semicon- 
ductor nanowires by measuring the weak antilocalization 
effectF™ 

The various possibilities of spin control in two- 
dimensional electron gases and planar wire structures 
opened up by the Rashba effect have inspired us to an- 
alyze theoretically the spin dynamics in tubular conduc- 
tors. We have used a cylindrical InAs nanowire with a 
surface two-dimensional electron gas as a model system, 
but our findings also apply to other systems, e.g. InN 
or InSb nanowires. In Sect. [II] we analyze the electronic 
states, focusing on spin properties, and we discuss the 
conditions under which a spin precession can be obse rved 
in tubular nanowires at zero ma gnet ic field (Sect. [TTT| ) and 
in an axial magnetic field (Sect.|TV|). In Sect.|V| we com- 
ment on the suitability of tubular conductors for spin 
electronic devices. 



II. ELECTRONS IN CYLINDRICAL WIRES 

Electrons confined in a cylinder move along the axis 
with a linear momentum hk (k real) and around the axis 
with an angular momentum HI (I integer). As long as the 
translational and rotational symmetries of the cylinder 
are not perturbed these momenta are conserved quanti- 
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ties. The wave function of an electron 

ip = exp (ikz) exp (il<fi) f (r) 

is a product of exponential functions in z, <fi, the coor- 
dinate along the axis and the azimuthal angle around 
the axis respectively, and a radial distribution function 
f (r). The distribution is determined by internal forces 
produced by the cylinder material. In our case, we took a 
planar 2-dimensional electr on gas (2DEG) at the surface 
of In As as a reference,* 31 * 32 * i.e. assuming a surface state 
charge density of Ns = 1.27 x 10 11 cm -2 , a background 
p-doping of rid = 2.8 x 10 17 cm -3 and an effective elec- 
tron mass of m* = 0.026 m e . The calculations were done 
for a cylinder radius ro = 50 nm. A schematic illustra- 
tion of the nanowire is depicted in Fig. [I] (upper inset). 
Electrons of atoms at the surface may find energetically 
more favorable states in the conduction band. Due to 
the Coulomb attraction between the electrons and the 
ions remaining at the surface the electrons get trapped 
in a layer close to the surface forming a 2DEG! 3 - 3 -^ The 
potential V resulting from the charge density of occu- 
pied electron states ipi :(T ,k, of ions at the surface and of 
dopants pBG 



l,a,k 



PBG /<* 
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is shown in Fig. [T] e is the elementary charge, a the spin 
index. e r — 14.6 is the bulk dielectric constant of InAs) 3 -^ 
It takes the polarization charges of the medium into ac- 
count. The potential profile is determined by Poisson's 
equation which is solved in cylindrical symmetry analyt- 
ically 



f r r' 
V= 47re e / r'dr'p (r') In — . 
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Equations ([!]) and ([2| are solved self-consistently. Start- 
ing from the potential of a homogeneous distribution of 
electrons in the cylinder the distribution is recalculated 
using the Schrodinger equation given below [see Eq. ([5|] 
and Eq. 0. The iteration procedure converges monoton- 
ically. We assumed an interface barrier of infinite height. 

Due to the electric field £ — —W/e across the surface 
of the cylinder the spin a of the electron is coupled to its 
orbital motion 
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The coupling-strength 7 is determined by the band struc- 
ture of the cylinder material (1.17 nm 2 for InAs)) 3 -^The 
second part of Eq. ([3| expresses %so in terms of Pauli 
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FIG. 1: Squared amplitude of the wave function 1^1 2 , the 
spinor components / and ft and potential profile V as a func- 
tion of the normalized radius r/ro. The upper inset shows a 
schematic illustration of the nanowires, including the relevant 
electric and magnetic fields. The lower inset shows the spin 
orientation along the circumference for j = 1/2. 



matrices for <J x ,y,z acting on a 2-component (spinor) wave 
function (^,^). The off-diagonal terms in Hso raise 
(lower) the value of the orbital angular momentum C z of 
^ (ipi) by ft. The stationary states are eigenstates of 
the total angular momentum J z — C z +S z (S z = ha z /2) 
with eigenvalues j = / zb 1/2. The spinor is of the form: 



e %kz e il<f> 



fir) 
te^h (r) 



(4) 



where /, ft are real functions and solve the differential 
equations 



2m 
ft 2 



k-yV ft 



- (h" + h'^j + - c) ft = k 7 v f . (5) 



Here, V h ± = (ft/) 2 /(2m*r 2 ) + V ± jV'l/r contains the 
contributions of the centrifugal force and the diagonal 
spin-orbit term, e = e — (ft/c) 2 /(2m*) is the energy with- 
out the axial kinetic energy. At the wire boundary we 
assumed a barrier of infinite height The influence of 
an external magnetic field B is not included yet. 

In Fig. [2] the energy e is plotted for several j-bands at 
B = 0. The parabola indicates the axial kinetic energy 
left out. It crosses the bands at the Fermi- momentum fcp, 
i.e. states with energy below the parabola are occupied. 
At k = the coupling between I and I + 1 vanishes [cf. 
Eq. ([3| ]. Classification with respect to I is possible. 
The splitting between the second and third band (I = 
±1) is caused by the diagonal part of %so and increases 
proportional to I for the higher states. Due to the mirror 
symmetry z — z states with angular momentum and 





FIG. 2: Energy vs. k dispersion at B = 0. The dashed line 
indicates the axial kinetic energy left out, which crosses the 
bands at k F - The pair of dots represent states forming the 
superpositions states ^-5/2, ± and ^+7/2, ± 5 while the square 
indicates the state ^+1/2, The dotted line illustrates the 
Gaussian wave packet of width 5k — 1/vq. 



FIG. 3: Spin density (st,s z ) of the lower energy states for 
total angular momenta j — 1/2,3/2,5/2 and 7/2. The spin 
is oriented only tangent ially and along the z-axis. 



spin reversed have the same energy. Therefore, all bands 
are twofold degenerate. 

The solution (/, h) of Eq. ^ for j = 1/2 at k F is shown 
in Fig. [l] The spin-orbit coupling increases linearly with 
fc, i.e. at k F with I = there is the strongest spin-orbit 
coupling. The spin density attains a sizable tangential 
component 

-(;i)(.^'r)(j;)=* 

The component along the wire axis is 

•-(8)0 axs)-'-*- ^ 

The radial component is zero. The spin orientation 
around the cylinder for j = 1/2 is illustrated in Fig. [I] 
(lower inset). According to Eq. (|6| the spin turns to the 
axial direction. This is shown for different values of j in 
the plot of the spin densities St and s z in Fig. [3j As can 
be seen here, the spin is oriented exclusively tangentially 
and along the axial direction. When averaged over the 
cylinder plane (• • • ) for each state ipj the spin compo- 
nents (a x ) and (a y ) are zero, while a finite contribution 
(<j z ) remains along the z-direction. 

III. SUPERPOSITION STATES AND SPIN 
PRECESSION 

For each k and j there are two solutions of Eq. ([5| 
ipj t ±. The (+)-state originates from (^,0), the (— )-state 



from (0,^+i), the solutions at k = 0, with j = I + 
1/2. They are orthogonal to each other and have opposite 
spin direction (±). They have different energies e and 
therefore different k F . Their superposition i/jj^ = (^ }+ + 
ijjj-)/\/2 yields 

(On = + 

+ (fj,+fj- - h 3,+ h 3~) cos ( k F,+ ~ k F -) z . 

The contributions of the basis states (±) almost cancel 
each other and are neglected further on. The interference 
between the states is constructive due to orthogonality 
and leads to 

(<7z)\\ ~ 2(/ j)+ / 7 -_)cos(fc F ,+ - k F -)z , (7) 

an oscillation of the average spin along the cylinder axis 
with a wavelength Ay = 27r/|fcp j + — &f,-|- The spin com- 
ponents (<7;e)||j (&y) || in the cylinder plane are both zero. 

Superpositions of eigenstates with different j's form 
states with a non-zero average spin component in the 
cylinder plane, e.g. i/jj^± = + ipj-i,-)/ V%- As 

one can easily retrace, these states originate from states 
with the same angular momentum I. The interference 
term gives the only ^-independent contribution to the 
densities of cr x ^a y . With 

( fir) \ 
ip ji+ = exp (izkF,+) exp {ilcj)) I ^ J , 

^-1,- = exp (izk F ^ exp [1 (I - 1) 0] ^ JJ(^ ) > 

(8) 
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the averages are 



(<7x)± = (hf) sin (k Fi + - fc F ; 
(cr 2/ )_ L = (ft/) cos (k Fi + — kp —^j 



(9) 



The (cr z )j_ contribution is small and does not depend on 
z. In particular, for the superposition ^1/2, ± of the lowest 
two states (<J Z }± is zero. 

For '0_5/2,_l, the superposition of ^_ 5 / 2 ,+ and ^_ 7 / 2 ,-, 
the spin precesses counterclockwise in the cylinder plane 
along the cylinder axis, as illustrated in Fig. |4ja). Here, 
we assumed an initial spin orientation along the — y di- 
rection, which in practice can be realized by spin injec- 
tion from a spin-polarized electrode. For ^+7/2, _l con - 
stituted of the opposite states ^+7/2,+ an d ^+5/2,- the 
spin precession is clockwise. Both precessions have the 
same period of X± = 27r/|/ci? + — fcp,-|- Their energy 
is degenerate. Due to their exactly inverse precession 
sense the combination of these states results in an oscil- 
latory behavior of the net spin orientation, as depicted in 
Fig. Qb). For an initial spin orientation along the —y di- 
rection the spin oscillates in the ?/z-plane. Superposition 
of the respective opposite states restores the left-right 
symmetry and eliminates spin precession. The oscilla- 
tion period A^ of ipj,± depends on j. For smaller e.g. 
'0-3/2, _l the corresponding difference in kp,+ and fcp 5 _ 
becomes smaller so that the period A^ is enlarged, as 
one can infer from Fig. [4jc) as compared to Fig. [4^b). 
The superposition state ^+1/2, ± constituted of the two 
lowest lying energy states ^±1/2, ± (cf. Fig. |2j square) 
shows no precession at all, because here fcp 5+ and kp,- 
are identical. Figure 4^d) shows the spin variation for 
a Gaussian wave packet of width 5k = 1/vq centered 
between the hp's of the states ^-3/2, ± and ^+5/2, _L- In 
position space this corresponds to a distribution of width 
2ro- The oscillation deviates from a purely harmonic os- 
cillation, as shown in Fig. |4jc), due to the contributions 
of the other states at the Fermi energy. This effect is also 
increasing with decreasing \j\ when the hp's get closer to 
each other. 

The electron spin is usually injected from a spin- 
polarized electrode in all states at the Fermi energy Ep 
having the correct spin direction. Thus, if only the di- 
rection of the spin is fixed by the electrode, all states are 
likely to transport electrons through the cylinder and a 
definite precession will not be observed. The total spin 
will only vary in the plane which is defined by the ini- 
tial spin orientation and the z-axis, similar to the situ- 
ation illustrated in Fig. [4j In order to observe spin pre- 
cession about the cylinder axis, a selection mechanism 
which breaks the left-right symmetry of the system must 
be adopted. As it will be discussed in the next section, 
this is achieved by applying a longitudinal magnetic field 
£=(0,0,fl). 




FIG. 4: (a) Counter-clockwise spin precession of electrons in 
the superposition state ^-5/2, ± at the Fermi energy consti- 
tuted of the states ^-5/2,+ and ^-7/2,- for a propagation 
along the wire axis from z/ro = to 30. (b) Spin orientation 
of the sum of the contribution shown in (a) and the corre- 
sponding clockwise contribution ^+7/2, _l being a superposi- 
tion of ^+7/2,+ and '0 +5 /2,-. (c) Spin oscillations resulting 
from the combinations of the two lower energy superposition 
states ips/2,± and ^+5/2, j_- (d) Spin variation for a Gaussian 
wave-packet of width 1/ro centered between the /cf's of the 
states ips/2,± and ^+5/2, _l (cf. Fig. [2J. 



IV. SPIN PRECESSION IN A MAGNETIC 
FIELD 



The vector potential A = (-By/2,Bx/2,Q) of a 
longitudinal magnetic field introduces a paramagnetic 
(Zeeman-) and diamagnetic (Landau-) term into Eq. ([5|. 
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FIG. 5: Energy vs. B dispersion (left panel) at k = and 
energy vs. k dispersion (right panel) at B = 0.13 T. The 
dashed line indicates the axial kinetic energy left out, which 
crosses the bands at Jgf- The pairs of dots indicate the states 
forming the superposition states ^-5/2, ± and ^+7/2, _l a t &f 
with a net spin in the cylinder plane. The two states ^+5/2,+ 
and ^+5/2,- with j — +5/2 are marked by triangles. 



V^± is extended to 

he n d ■ gm*\ e 2 B 2 2 

Vi ± = Vi ± + B I ± + r 2 , 

2m* \ 2m e J 8m* 

with g the gyromagnetic-factor of the electron spin 
( — 14.9 for InAa^l). The paramagnetic (second) term 
in Vi y ± raises e for states with j(or I) > and lowers 
e for states with j(or I) < 0. The energy difference in- 
creases ex IB for B <C ifi/(ero) (cf. Fig. |5j. For larger 
B e increase ex B 2 due to the diamagnetic (third) term. 
In the linear range the influence of the r-dependence of 
the third term is negligible. The densities do not change 
significantly. 

The main effect of B is the energetic separation of 
the ±j-states. It opens possibilities of observing spin 
dynamics in electronic transport. This will be demon- 
strated in the following at B = 0.13 T. Figure [5] shows 
the independence at k = up to B = 0.13 T and the k- 
dependence at B = 0.13 T of e for states from j = ±1/2 
to ±9/2. Again, the parabola marks the Fermi edge. Su- 
perpositions with spin in the cylinder plane according to 
Eq. ([£]), i/jj : ± are marked as pairs in Fig. [5] The lower pair 
corresponds to ^-5/2, ± depicted in Fig. |4^a) . As illus- 
trated in Fig. [6ja), it shows the same counter-clockwise 
precession. In contrast to the zero field case, now the su- 
perposition state ^+7/2, ± has a larger kp -difference, i.e. 
a shorter precession length [cf. Fig.[6jb)]. Consequently, 
the precessions of ^-5/2, _l an d ^+7/2, _l are n °t exactly 
opposite. In contrast to the case at B — 0, the spin 
still rotates following the state with the faster precession, 
when both states are superposed. This is illustrated in 
Fig. [6^c), where one finds that in addition to the oscilla- 
tion of the spin amplitude its orientation is also changed 
during propagation. Thus, by applying a magnetic field 




FIG. 6: (a) Spin precession of electrons at the Fermi energy 
propagating along the wire axis for the superposition state 
'0-5/2 > _l at B = 0.13 T. (b) Corresponding spin precession 
for the state ^+7/2, ±- (c) Spin orientation and magnitude 
of the sum of the contributions shown in (a) and (b) for a 
propagation from z/tq — to 30. The arrow indicates the 
direction of the initially injected spin. 



a spin precession can be achieved. 

In the previous section, we already pointed out that 
the superposition state ipju with equal total angular mo- 
mentum but opposite spin orientation result in an os- 
cillation of the average spin along the cylinder axis. In 
Fig. FF|a) and (b) these oscillations of {cr z )\\ are shown 
for ainerent values of j at B = 0.13 T. One finds that 
for larger total angular momentum values the oscillation 
period is shorter owing to the larger difference of Fermi 
wave vectors. In Fig. [5] the states contributing to ^+5/2,11 
are marked by triangles. Compared to the previously dis- 
cussed i/jj^ states, here the difference in the Fermi vec- 
tors is relatively large, leading to a faster oscillation com- 
pared to the spin precession period shown in Fig. |6jb). 

Once again the application of an axial magnetic field 
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z/r j=+5/2 



FIG. 7: (a) Spin orientation (a z ) along the wire axis for the 
superposition states ipj^ with j = +1/2, +3/2 and +5/2 at 
B — 0.13 T. (b) Illustration of the oscillation of the average 
spin along the wire axis for the superposition state ^+5/2,11- 
(c) Modulation of the total spin orientation (a z ) (dashed line) 

resulting from a combination of the *0 5/2, || an d ^+5/2, || states 

at B = 0.13 T. 



breaks the symmetry of the ijj±j\\ states. As can be 
inferred from Fig. (7^c), a different oscillation period is 
found for the ^+5/2, || an d '0-5/2, || states. Thus, when 
these states are combined a beating in the oscillation of 
the average spin appears. 

V. CONCLUSIONS 

In the previous two sections we learned that an injected 
spin is strongly modulated while propagating through 
a cylindrical nanowire. For a spin injection along the 
wire axis, e.g. by a ferromagnetic electrode, the spin 
is carried by superposition states with equal total angu- 
lar momenta. In analogy to the spin field-effect tran- 
sistor based on a planar 2DEG, 10 a transistor structure 
can be realized by placing a second magnetic electrode 
at the opposite terminal of the nanowire as a spin de- 
tector. Control of the spin orientation can achieved by 
manipulating the strength of the Rashba effect by means 
of a gate electrode. By applying a bias voltage to the 
gate, the strength of the electric field £ in the surface 
2DEG is adjusted. In order to obtain a uniform control 
within the channel, a so-called wrap-around gate should 
be preferred] Usually, in a realistic situation a larger 
number of states with different total angular momenta 



j is occupied. As we observed, for each superposition 
state 7/^,11 different oscillation periods are found. This 
leads to a rather complex modulation of the spin along 
the axial direction. An obvious strategy for simplifica- 
tion is to reduce the number of occupied states, i.e. by 
depleting the channel by means of a gate. Another pos- 
sibility might be to only occupy certain states by means 
of /c-selective filters. This might be realized by means 
of an injection through a single or a resonant tunneling 
barrier. As pointed out in Sect. |III[ one possible way 
to model this situation is to assume the formation of a 
state with a Gaussian distribution around the average 
momentum. 

In addition to a spin injection and detection along the 
wire axis it is also possible to inject spins in transversal 
direction. Here, the spins are carried by superposition 
states ipj,±_ constituted of states with different total an- 
gular momenta j. As long as no magnetic field is applied 
the spin is exclusively modulated in the plane spanned 
by the injection orientation and the wire axis. Here, the 
output signal in a spin field-effect transistor is gained by 
gate-modulating the spin orientation along or opposite 
to a detector electrode, which is polarized parallel or an- 
tiparallel to the injector. By applying an axially oriented 
magnetic field, spin precession about the wire axis can be 
achieved. This additional feature might be an interesting 
option to implement more complex functionalities in spin 
electronic devices. 

In conclusion, we have shown that semiconductor 
nanowires affected by Rashba spin-orbit coupling are 
promising candidates for future nanowire-based spin elec- 
tronic devices. The complex spin dynamics in these 
cylindrically-shaped conductors provide many opportu- 
nities to tailor the device functionality. 
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